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Abstract 

We present results from a 150 ksec Suzaku observation of the Seyfert 1 NGC 
3516 in October 2005. The source was in a relatively highly absorbed state. Our 
best-fit model is consistent with partial covering by a lowly-ionized absorber with a 
column density near 5xl0 22 cm -2 and with a covering fraction 96-100%. Narrow 
K-shell absorption features due to He- and H-like Fe confirm the presence of a high- 
ionization absorbing component as well. A broad Fe Ka diskline is required in all fits, 
even after the complex absorption is taken into account; an additional partial-covering 
component is an inadequate substitute for the continuum curvature associated with 
the broad line. The narrow Fe Ka line at 6.4 keV is resolved, yielding a velocity width 
commensurate with the optical Broad Line Region. The strength of the Compton 
reflection hump suggests a contribution mainly from the broad Fe line origin. We 
include in our model soft band emission lines from He- and H-like ions and radiative 
recombination lines, consistent with photo-ionization, though a small contribution 
from collisional ionization is possible. 

Key words: galaxies: active — X-rays: galaxies — galaxies: individual (NGC 
3516) 


1. Introduction 

In Seyfert AGN, the 6.4 keV Fe K a emission line has long been known to be an important 

diagnostic of the material accreting onto supermassive black holes. The associated Compton 

reflection hump, frequently seen in Seyfert spectra above ~7 keV and peaking near 20-30 

keV (Pounds et al. 1990), indicate that Seyferts’ Fe lines may have an origin in optically 

thick material, such as the accretion disk. Observations with ASCA indicated that many 

Fe Ka lines were broad (FWHM velocities up to ~0.3c) and asymmetrically skewed towards 
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lower energies, implying an origin in the inner accretion disk; the line profile is sculpted by 
gravitational redshifting and relativistic Doppler effects (e.g., Tanaka et al. 1995, Fabian et al. 
2002). However, XMM-Newton and Chandra observations have been revealing a more complex 
picture. A narrow Fe K component (FWHM velocities ~5000 km s -1 or less) instead appears 
to be much more common; these lines’ widths suggest emission from distant material, such 
as the outer accretion disk, the optical Broad Line Region or the molecular torus. Spectral 
observations in which the broad and narrow components are deconvolved are thus a prerequisite 
for using the Fe K line as a tracer of the geometry of the emitting gas. 

At the same time, there is strong evidence from X-ray and UV grating observations for 
the presence of ionized material in the inner regions of a large fraction of AGN (e.g., Blustin et 
al. 2005). High-resolution spectroscopy shows the gas is usually outflowing from the nucleus; 
typical velocities are ~ a few hundred km s -1 . Absorption due to a broad range of ionic species 
is commonly seen; and for many sources, there is evidence for several different photoionized 
■absorbing components, as opposed to a single absorber, along the line of sight. In the Fe K 
bandpass, some Seyferts show evidence for absorption by H- or He-like Fe, indicating a zone of 
highly-ionized absorbing material (e.g., NGC 3783, Reeves et al. 2004). 

Cold absorbing gas, with line of sight columns in excess of the Galactic value, routinely 
observed in Seyfert 2 AGN in accordance with unification schemes (Urry, Padovani 1995), 
have also been reported in some Seyfert 1 AGN. Importantly, variations in column density 
on timescales from hours to years have been observed in both Seyfert 1 AGN (e.g., Lamer et 
al. 2003, Puccetti et al. 2004) and Seyfert 2 AGN (Risaliti et al. 2002; Risaliti et al. 2005), 
suggesting that the absorbing circumnuclear material is not homogeneous in either Seyfert type. 

NGC 3516 is a well-studied, nearby (z=0. 008836; Keel 1996) Seyfert 1 AGN that can 
display strong 2-10 keV flux variability on timescales of hours to years (e.g., Markowitz, Edelson 
2004). Previous X-ray spectroscopic observations of NGC 3516 e.g., Nandra et al. (1997) using 
ASCA , have indicated the presence of a broad Fe K line, but this source is known to also 
contain complex and ionized absorption. Numerous UV absorption lines, including N V, C IV 
and Si IV, were observed with the International Ultraviolet Explorer (Ulrich, Boisson 1983); 
absorption line strengths vary on timescales as short as weeks as the absorber responds to 
variations in the ionizing flux (e.g., Voit et al. 1987). Hubble Space Telescope observations 
have revealed that this component of absorbing gas (henceforth called the “UV absorber” ) may 
consist of several distinct kinematic components (Crenshaw et al. 1998). 

X-ray spectra of NGC 3516 can exhibit evidence for large columns (>10 22 cm -2 ) of 
absorbing gas (e.g., Kolman et al. 1993), though the X-ray absorbers can also display variability 
on timescales of years (e.g., Mathur et al. 1997). Using Chandra gratings data from observations 
in April 2001 and November 2001, Turner et al. (2005) observed K-shell absorption lines due 
to H- like Mg, Si and S, and He-like Si, evidence for a highly-ionized absorber, likely with 
column density >10 22 cm -2 , outflowing at ~1100 km s -1 . Simultaneous with these Chandra 
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observations in 2001 were two XMM-Newton observations. Turner et al. (2005) modeled the 
continuum curvature of the two XMM-Newton EPIC spectra by including a partial covering, 
mildly-ionized absorber; the column density was ~2.5xl0 23 cm -2 , with a covering fraction of 
~50%. However, a direct consequence of including this component in the spectral model was 
that the formal requirement for the broad Fe line vanished, leading to uncertainty as to whether 
the broad Fe line really existed in NGC 3516. Spectral fitting using an instrument with a wide 
bandpass is thus necessary to remove such model degeneracies. 

In this paper, we report on an observation of NGC 3516 made with the Suzaku obser- 
vatory in October 2005. The combination of the X-ray Imaging Spectrometer (XIS) CCD and 
the Hard X-ray Detector (HXD) instru m ents have yielded a broadband spectrum covering 0.3 
to 76 keV, allowing us to deconvolve the various broadband emitting and absorbing compo- 
nents. Furthermore, the exceptional response of the XIS CCD and high signal-to-noise ratio 
of this observation have allowed us to study narrow emission lines in great detail. §2 gives a 
brief overview of the Suzaku observatory, and describes the observation and data reduction. 
§3 describes fits to the time-averaged spectrum. Variability analysis is briefly discussed in §4. 
Flux- and time-resolved spectral fits are discussed in §5. The results are discussed in §6, and a 
brief summary is given in §7. 

2. Observations and data reduction 

The nucleus of NGC 3516 was observed by Suzaku from 2005 October 12 at 13:45 UT 
until October 15 at 09:07 UT. Suzaku was launched 2005 July 10 into a low-Earth orbit. It 
has four X-ray telescopes (XRTs; Serlemitsos et al. 2007), each with a spatial resolution of 2' 
(HPD). The XRTs focus X-rays onto four X-ray Imaging Spectrometer (XIS; Koyama et al. 
2007) CCDs, which are sensitive to 0.2-12 keV X-rays on a 18' by 18' field of view, contain 
1024 by 1024 pixel rows each, and feature an energy resolution of ^150 eV at 6 keV. Three 
. CCDs (XIS0, 2 and 3) are front-illuminated (FI), the fourth (XISl) is back-illuminated (BI) 
and features an enhanced soft X-ray response. The XRT /XIS combination yields effective areas 
per detector of roughly 330 cm 2 (FI) or 370 cm 2 (BI) at 1.5 keV, and 160 cm 2 (FI) or 110 cm 2 
(BI) at 8 keV. Each XIS is equipped with two 50 Fe calibration sources which produce fluorescent 
Mn Ka and K/3 lines and are located at the CCD corners. Suzaku also features a non-imaging, 
collimated Hard X-ray Detector (HXD; Takahashi et al. 2007); its two detectors, PIN and GSO, 
combine to yield sensitivity from ~10 to ~700 keV. Further details of the Suzaku observatory 
are given in Mitsuda et al. (2007). 

2.1. XIS Reduction 

The XIS data used in this paper were version 1.2 of the screened data (Fujimoto et al. 
2007) provided by the Suzaku team. The screening is based on the following criteria: grade 
0, 2, 3, 4, and 6 events were used, the cleansis script was used to remove hot or flickering 
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pixels, data collected within 256 s of passage through the South Atlantic Anomaly (SAA) were 
discarded, and data were selected to be 5° in elevation above the Earth rim (20° above the 
day-Earth rim). The XIS-FI CCDs were in 3x3 and 5x5 editmodes, for a net exposure time 
after screening of 135.0 (XISO), 134.8 (XIS2) and 135.2 (XIS3) ks. XIS1 was also in 3x3 and 
5x5 editmodes, for a net exposure of 135.4 ks. The XIS was in normal clocking mode. 

The source was observed at the nominal center position of the XIS. For each XIS, we 
extracted a 3' radius centered on the source. The background was extracted using four circles 
of radius 1.5', each located ~6' from the source, but chosen to avoid the z=2.1 QSO RX 
J 110741.4+723235, located 4.5' SE of NGC 3516. Spectra were binned to a minimum of 50 
counts bin -1 to allow use of the x 2 statistic. 

Response matrices and ancillary response files (ARFs) were generated for each XIS in- 
dependently using XISSIMRMFGEN and XISSIMARFGEN version 2006-10-26 (Ishisaki et al. 2007). 
The ARF generator takes into account the level of hydrocarbon contamination on the optical 
blocking filter. We estimate a carbon column density of 0.8, 1.2, 1.7, and 2.8 xlO 18 cm" 2 for 
XISO, 1, 2 and 3, respectively. Finally, we co-added the three XIS-FI spectra using mathpha, 
and co-added the response files and ARFs using addrmf and ADDARF, respectively. 

To examine the accuracy of the XIS RMFs and determine residual line width due e.g., 
to imperfect CTI correction, we generated spectra for the 55 Fe calibration source lines on each 
XIS using the above response matrices and ARFs. We fit the calibration source spectra with 
three Gaussians. Two Gaussians were for the Mn Ka doublet (expected energies 5.899 keV 
and 5.888 keV), with energy centroids fixed to be 11 eV apart, and the higher energy line flux 
set to twice that of the lower energy one. The third Gaussian was used to model the K(5 line, 
expected at 6.490 keV. We found the average of all the calibration line widths a to be < 4 
eV. The Mn Ka line energy centroids for the co-added FI spectrum were consistent with the 
expected energies to within 1 eV. For XISl, the Mn Ka line energy centroids were 3±2 eV 
lower than expected. Such discrepancies are well within the accuracy (~0.2% at Mn-Kcc) of 
the energy calibration of XIS. Fitting the calibration source lines without the response file, we 
determined the FWHM energy resolution during the observation to be 145 eV (average of the 
4 XlSes). 

2.2. HXD Reduction 

We used data from the HXD-PIN only; NGC 3516 was extremely faint in the HXD-GSO 
band, and some aspects of the GSO background are still being studied, so we defer analysis of 
the GSO data to a later time. The PIN source spectra were extracted from cleaned version 1.2 
(pre-1.2-rl) HXD event files provided by the HXD instrument team. PIN background count 
rates are variable and strongly depend on the time since SAA passage (Kokubun et al. 2007), so 
we selected data according to the following criteria: at least 500 s since SAA passage, COR > 8 
GV, and day- and night-Earth elevation angles each >5°. Instrumental (non-X-ray) background 
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spectra for the PIN were provided by the HXD Team (“Background A” model) generated from 
a t im e-dependent model. The model utilized the count rate of upper discriminators as the 
measure of cosmic ray flux that passed through the silicon PIN diode and yielded background 
spectra based on a database of non X-ray background observations with the PIN (Fukazawa 
et al. 2007). The current accuracy of the PIN NXB model for a 1 day observation is about 
5 % (peak-to-peak residuals). Both the source and background spectra were generated with 
identical good time intervals, and the source exposure was corrected for instrument dead time 
(a ~5% effect). This yielded a good time exposure of 105.8 ks. 

Data <12 keV were discarded due to noise contamination near the lower threshold of 
the PIN diode. Data above 76 keV were also discarded: the gain above an internal Bi Ka 
calibration line at 76 keV is not well-defined, though the photon statistics above this energy 
were poor anyway for this observation. Further details of the HXD in-orbit performance are 
given in Kokubun et al. (2007): To model the contribution to the total background from 
the Cosmic X-ray Background (CXB), a spectrum of the form 9.0xl0" 9 (E/3keV) _o ' 29 exp(- 
E/40keV) erg cm -2 s -1 sr -1 keV -1 (Gruber et al. 1999) was used. We note that some recent 
works (e.g., Prontera et al. 2007) suggest a 10% normalization increase compared to Gruber et 
al. (1999). However, the effect on the net spectrum was negligible; for instance, the change in 
Compton reflection component strength was 1%, typically. To simulate a CXB spectrum using 
XSPEC, we assumed a model of the form powerlaw*highecut, with photon index r=1.29, 
and low- and high-energy cutoffs of 0.1 and 40 keV, respectively. The power-law normalization 
of 8.8 xlO -4 ph cm -2 s" 1 keV -1 (at 1 keV) was used, appropriate for a source observed in 
XlS-nominal mode. The 12-76 keV CXB flux was 1.4xl0~ u erg cm -2 s _1 (using the Gruber 
et al. 1999 normalization). The total (X-ray plus particle) background 12-76 keV flux was 
4.4xl0~ 10 erg cm -2 s _1 . 

The spectrum was binned to a minimum of 400 count bin -1 . We used the response 
file aeJixd_pinxinom_20060814.rsp. The mean 12-76 keV net source flux and count rate were 
l.lxlO -10 erg cm -2 s -1 and 0.16 ct s -1 . Figure 1 shows the net source, background, and total 
(source + background) spectra. The source spectrum is always at least 15% of the total up to 
~40 keV. 

3. Model Fits to the Time- Averaged Spectrum 

We used 0.4-1. 5 keV data in the XIS-FI spectrum and 0.3-10.0 keV data in the XIS-BI 
spectrum. We ignored 1.72-1.87 keV in the co-added FI spectrum and 1.80-1.87 keV in the 
BI spectrum due to uncertainties in calibration associated with the instrumental Si K edge. 
In all fits, we included a constant to account for relative instrument normalizations. We left 
the relative XIS-BI/XIS-FI normalization free, but best-fit values were always extremely close 
to 1.00. The PIN/XIS-FI normalization was kept fixed at 1.13, a value derived using Suzaku 
observations of the Crab (Yamasaki et al, in prep.). All errors on one interesting parameter 
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correspond to Ay 2 = 2.71 (with the XIS-BI/XIS-FI normalization left free). A redshift of 
z=0. 008836 was used throughout. The abundances of Lodders (2003) were used. A neutral 
Galactic column of 2.94xl0 20 cm -2 was included (Dickey, Lockman 1990). 

3.1. Preliminary Broadband Fits 

The X-ray continuum emission of Seyfert 1 AGN is usually dominated by a power-law 
component thought to originate from inverse Comptonization of soft seed photons by a hot 
corona near the central black hole (e.g., Shapiro et al. 1976; Sunyaev, Titarchuk 1980; Haardt 
et al. 1994). A simple power-law (henceforth denoted the “primary power-law”) over 0.3-76 keV 
yielded a very poor fit, with x 2 /do/ (degrees of freedom) = 49708/1449. As shown in Figure 
2(a), residuals strongly indicated the need to include absorption to account for continuum 
curvature <3 keV. The narrow 6.4 keV Fe Ka line was also obvious. 

We therefore added an absorbing column to the power-law, assuming a covering fraction 
of 1 and systemic redshift (this component is henceforth denoted the “primary absorber”). We 
used the xstar table grid25BIG-mt.fits. We also added Gaussians to model Fe Ka and K (3 
emission lines. The Fe K/3 energy centroid was kept fixed at 7.06 keV; the normalization was 
kept fixed at 0.13 times that of the Ka line. With these changes to the model, y 2 /do/ fell to 
8311.962/1444. The best-fit ionization parameter and column density were log(£) = 2.0±0.1 
erg cm s _1 and 7.1A0.2 xlO 22 cm -2 , respectively (£ = 47rFi on /n; F ion is the 0.5-20 keV ionizing 
continuum flux; n is the density of the reflecting material). However, as shown in Figure 2(b), 
this model did not accurately describe the broadband emission. 

To attempt to model the soft emission, we added a second power-law (the “soft power- 
law”), with photon index V tied to the primary power-law; ^fdof fell to 3038.86/1443. As 
shown in Figure 2(c), residuals suggested the presence of soft X-ray emission lines, e.g., near 
0.56 keV, likely due to O VII. 

We added 13 Gaussians to our model, widths a were fixed at 0.5 eV. Energy centroids 
for the lower signal-to-noise ratio lines were kept fixed at lab-frame energies. Table 1 shows the 
results for the lines in our best-fit baseline model (see below). Data/model residuals for these 
emission lines are shown in Figure 3. It was significant at > 99.0% confidence in an F-test to 
include each line in the fit. We identify these lines as originating in H-like C, N, O, Ne and 
Mg, and He-like N, O, Ne and Mg. We also report radiative recombination (RRC) lines of 
O VII, O VIII and Ne IX (and possibly C VI, blended with the 0.500 keV N VII emission 
line). The lines are likely attributed to photo-ionization. Most of these lines have been reported 
previously; we refer the reader to Turner et al. (2003) for results using the XMM-Newton- RGS. 
The strongest line detected in both the RGS spectrum and the Suzaku spectrum is O VII. We 
also included a line near 0.83 keV for Fe L XVII, to model any contribution from collisionally- 
ionized gas. Turner et al. (2003) included a Mg XI recombination edge component; we eschew 
it in favor of a Mg XII line at 1.47 keV. Including these emission lines in the model resulted 
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in x 2 /dof falling to 2346.72/1424. Figure 2(d) shows the residuals to the data after the soft 
X-ray lines are modeled. 

Residuals in the PIN band, peaking near 20-30 keV, signaled the need to include a 
Compton reflection component. We added a pexrav component, assuming solar abundances, 
a cut-off energy fixed at 400 keV, an inclination fixed at 30°, and a normalization tied to 
that of the primary power-law. Initially, we assumed that the reflector is subject to the same 
absorption as the primary power-law. In the best fit-model, x 2 /dof fell to 1895.49/1423, and 
the value of the reflection fraction R (defined as Q/ 2zr, where Q, is the solid angle subtended by 
the reflector) was 2.8lo;|. The photon index P was 2.142j]o;oo7- As shown in Figure 2(e), the 
residuals in the PIN band were thus corrected. 

The uncertainty on the relative PIN/XIS cross-normalization is about 1%; the subse- 
quent effect on R is smaller than that associated with the uncertainty of the PIN NXB. We 
modified the intensity of the PIN NXB by ±2%, which is the la level of the current repro- 
ducibility of the PIN NXB. R changed by ±0.2, which is smaller than the statistical error on 
R. Errors quoted on R for the remainder of this paper are the statistical errors only; readers 
should bear in mind the additional systematic uncertainty associated with the NXB. 

Finally, we re-fit the model, assuming that the PEXRAV component was not affected 
by the primary absorber. This yielded a goodness of fit nearly identical to the previous fit, 
with x 2 /dof — 1479.3/1417. All fit parameters were virtually identical to the previous fit; for 
instance, R was 1.9±0.4. 

There remained residuals in the 5-6 keV band. We first discuss modeling these residuals 
using relativistic diskline components, and later (§3.2) we will discuss if a partial covering 
component could provide as good a fit. We added two LAOR components (Laor 1991) for 
relativistic Fe Ka and Fe K j3 emission. Formally, the K/3 diskline is not required (omission 
of this component does not change x 2 /do f significantly), but we include it for completeness. 
The normalization of the K/3 diskline was fixed at 0.13 times that of the Ko- diskline. The Fe 
Ka line energy was constrained to lie within rest- frame energies of 6.40 and 6.96 keV; the K/3 
line energy was fixed at 7.06 keV. All other parameters were kept tied between the Ka and 
K/3 components. The emissivity index /3 (used when quantifying the radial emissivity per unit 
area as a power-law, r -/3 ) was fixed at 3.0. Initially, the disk inclination i was fixed at 30°. 
The outer radius R out was kept fixed at 400 R g (1 R g = GMbh/c 2 ). The inner radius i? in was 
left free but set at 2 R g initially. With the best-fit model, x 2 /dof fell to 1642.51/1420, and 
the data/model residuals near 5-6 keV were improved considerably; see Figure 2(f). It was 
significant at >99.99% confidence in an E-test to thaw z; in this and all subsequent models, we 
set the inclination of the pexrav component to that of the diskline. x 2 /dof fell to 1521.46/1419 
for i < 23°. We refer to this model as the “1-absorber + Compton reflection + diskline” model. 
Other model parameters are listed in Table 2. 

Visual inspection of the data/model residuals in the Fe K bandpass suggested some dips 
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near 6.7 and 6.9 keV, at the rest-frame energies for Fe XXV and Fe XXVI; see Figure 4(a). 
Such absorption features might be expected, given the detection of a high-ionization absorber 
by Turner et al. (2005). Adding a narrow (a = 0.5 eV) Gaussian with energy centroid fixed at 
6.70 keV yielded an improvement in fit of Ax 2 =-6.3 (for one less dof ), significant at 98.4% in 
an F-test. The absolute values of the intensity and equivalent width EW were 4±2 xlO -6 ph 
cm -2 s _1 and 9±4 eV. Adding a Gaussian at 6.96 keV, x 2 fell by 7.9, significant at 99.4% in 
an F-test. The absolute values of the intensity and equivalent width EW were 3±2 xlO -6 ph 
cm -2 s -1 and 9±5 eV. 

The Chandra-HETGS spectrum yielded narrow absorption features ~l-3 keV due to 
highly-ionized Mg, Si, and S, but we do not significantly detect any narrow absorption features 
at those energies. A small dip in the spectrum near 2.3 keV is close to the expected energy for 
S XII. However, XSTAR models demonstrate that S XV absorption, though not significantly 
detected here, is always stronger than S XII. This feature is more likely due to calibration 
uncertainty associated with an instrumental Au M edge. 

To model these absorption features, we added a second XSTAR component (henceforth 
denoted the “high-ionization absorber”). Based on the results of Turner et al. (2005), we 
assumed an outflow velocity of 1100 km s -1 . X 2 /do f fell to 1485.04/1417. In this model, 
the high-ionization absorber had a column 4.0jlf;f x 10 22 cm -2 , similar to the value used by 
Turner et al. (2005), and log(£) = 3. 71^7 erg cm s -1 . We henceforth refer to this model as 
our “2- absorber” or “baseline” model. Data/model residuals are shown in Figures 2(g) and 
4(b). The primary power-law, with T = lVlj^osj was absorbed by a column 5.5±0.2 x 10 22 
cm -2 and log(£)=0.3±0.1 erg cm s _1 . Forcing the ionization parameter to a much lower value 
(e.g., log(0 0.5) resulted in a significantly worse fit, with large residuals at 1. 5-3.0 keV; 
this is likely because the lower ionization does not account for absorption edges due to higher 
ionization species of Si and S. We caution, however, that this result could be influenced by 
residual calibration effects in the XIS near 1.8 keV and 2.3 keV (instrumental Si K and Au M 
edges, respectively). Thawing the diskline emissivity index (3 yielded no improvement to the 
fit. We found (3 was 3.2±2.9; we leave (3 fixed at 3.0 in subsequent fits. Figure 4(c) shows the 
data/model residuals when all narrow and broad lines are removed from the fit. Figure 4(d) 
shows the residuals to the best-fit baseline model when the diskline components are removed 
and the model re-fit; large residuals and a significantly higher value of x/dof (1650.66/1421) 
compared to the baseline model indicated that removing the disklines yielded a significantly 
worse fit (>99.999% in an F-test). In this model, the strength of the Compton reflection hump 
was R = 1.7lo;|; a contour plot of R as a function of T is shown in Figure 5. Also shown 
in Figure 5 is a contour plot of R as a function of the PIN /XIS-FI normalization, which had 
been kept fixed at 1.13; in this plot, F was a free parameter. Figure 6 shows an unfolded model 
spectrum. Other model parameters are listed in Table 2 (see Table 1 for the soft X-ray emission 
lines). 
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3.2. Additional partial covering components 

We next explored the possibility of an additional, partial covering component. Starting 
with the baseline model, we added a partial covering component consisting of a power-law 
(with photon index tied to that of the primary power-law) absorbed by an XSTAR component. 
We first kept the column density 1Vh,pc and ionization parameter log(£pc) tied to those of the 
primary absorber. In the best-fit model (henceforth denoted Model PCI), the added power-law 
had a normalization 0.32±0.03 times that of the primary power-law. As shown in Table 3, N h 
and log(£) of the primary absorber and the diskline parameters were consistent with values 
of the baseline model. However, x 2 /dof was 1485.69/1416, virtually identical to the baseline 
model, and so there is no formal requirement to include the new partial covering component 
when the column density and ionization states are tied to those for the primary absorber. 

Next, we untied Ah,pc and log(£pc) and refit (denoted Model PC2). The best-fit pa- 
rameter values are listed in Table 3. Compared to Model PCI, x 2 dropped by only 5.2 for 
2 additional degrees of freedom, significant at only ~90% in an F-test. It is thus not highly 
significant to include the new partial covering component in this case either. 

Next, we addressed whether it was possible for a partial covering component to mimic 
the curvature in the ~4-6 keV continuum modeled above as a relativistically-broadened Fe 
line. Starting with the baseline model, we removed the disklines, and added a partial covering 
component consisting of a power-law (with photon index tied to that of the primary power- 
law) absorbed by an XSTAR component with a relatively low value of the ionization parameter. 
The column density needed to be >10 22 cm -2 in order to have a cutoff near ~4-6 keV. In the 
best-fit model, the partial coverer had a column of 7.8^f 8 8 1 x 10 22 cm -2 , and log(£) = -0.4 
(errors pegged at both hard limits). The new partially-covered power-law had a normalization 
<0.22 times that of the primary power- law. However, x 2 /dof was 1676.57/1418, a much worse 
fit compared to the baseline model, and there were still large data/model residuals in the Fe 
K bandpass; see Figure 2(h). We conclude that a partial coverer cannot mimic the observed 
curvature of the diskline. 

3.3. Relativistic reflection fits 

Since the broad diskline is a signature of Compton reflection off a possibly-ionized disk, 
the ionized reflection models of Ross, Fabian (2005) can be used to model the entire reflection 
spectrum (broad Fe line plus reflection continuum) in a self-consistent manner. We removed 
the diskline and PEXRAV components from the baseline model, and instead used the grid of 
models of Ross, Fabian (2005); the ionization parameter £ was a free parameter but the photon 
index of the illuminating continuum was tied to that of the primary power-law. The entire 
reflection spectrum was then convolved with a relativistic LAOR profile; free parameters here 
were the emissivity index /?, inner disk radius R m and disk inclination i. We also included a 

1 parameter error pegged at lower limit 
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second, unblurred PEXRAV component to model any Compton reflection associated with the 
narrow 6.4 keV line. 

However, due to the heavy absorption present in NGC 3516, any putative soft X-ray 
excess resulting from blurred and ionized reflection would not be observed. The best-fit model, 
which had y 2 /do f = list chi-squared suggested a fairly neutral or lowly-ionized reflector, with 
log(0 < 1.8 erg cm s _1 . Parameters for the blurring disk were similar to what was obtained in 
the baseline model: an inclination of 25±10°, an inner disk radius < 4 R g , and an emissivity 
index of 3.0±0.5. The total Compton reflection strength R was 2.0 list uncertainty on R, with 
2/3 of that associated with the disk reflection, and 1/3 with the narrow Fe line. 

3-4- Narrow Fe K line properties 

In our baseline model, the best-fit energy centroid for the narrow Fe Ka line was 
6.398±0.004 keV, consistent with neutral Fe. The observed line width <r 0 bs was 261*3 eV. 
The intrinsic line width <7 intr was found by subtracting in quadrature the 55 Fe calibration line 
width a of <4 eV from the measured line width. We inferred a[ ntr = 261*4 eV, which cor- 
responds to a FWHM velocity width of 2800l*5oo km s~*. This width is consistent with the 
results obtained by CTmndra-HETGS for the two observations in 2001, 1290li29o an d 36301^4° 
km s -1 (Yaqoob, Padmanabhan 2004). 

Figure 7 shows a contour plot of the broad line intensity versus the narrow- line intensity, 
illustrating that the two lines are detected independently at >4cr confidence. Such a result is 
a product of the combination of the narrow response of the XIS (yielding extremely high sig- 
nal/noise in the narrow line) and Suzaku’s broad bandpass. Similar results have been reported 
e.g., for the Suzaku observation of NGC 2992 (Yaqoob et al. 2007). 

Finally, we discuss limits to a Fe Kcc Compton shoulder. We added a Gaussian emission 
line at 6.24 keV (rest-frame), with width tied to that of the Ka core. We found an upper limit 
to the intensity of 7xl0 -6 ph cm -2 s -1 , or 13% of the Ka core intensity. This limit corresponds 
to an EW of 21 eV. 

4. Timing analysis 

To compare the variability properties of the primary and soft power-laws, we extracted 
light curves, summed over all four XlSes, orbitally-binned, and background-subtracted, for the 
0. 3-1.0 and 2-10 keV bands. These light curves are plotted in Figure 8, along with the 12 — 76 
keV orbitally-binned and background-subtracted PIN light curve. For each light curve, we 
calculated the fractional variability amplitude F V3X (which quantifies the variability in excess 
of measurement noise) and its uncertainty following Vaughan et al. (2003). For the 2-10 keV 
band, F var was 9.3 ± 0.2%, with a maximum/minimum flux ratio of roughly 1.4. F var for the 
0. 3-1.0 keV was 1.3 ±1.8%. For the 12-76 keV band, no significant variability in excess of that 
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due to measurement errors was detected (F var undefined 2 ). We note that the systematic error 
associated with the non-X-ray background for a 5760 s bin was roughly 6 % for the 15-40 keV 
band (Mizuno et al. 2006). 

5. Flux-resolved spectral fits 

We performed flux-resolved spectral fits to search for any physical connection between 
the soft and primary power-laws. Despite the limited flux range exhibited during this obser- 
vation, we attempted to determine, e.g., if the observed X-ray flux variability could be due to 
rapid variations in column density of the primary absorber, if both power-laws vary together, or 
if one power-law is constant. We split the time-averaged spectrum into periods when the 2-10 
keV flux was higher and lower than the average 2-10 keV flux of 2.31 xlO -11 erg cm -2 s _1 , as 
illustrated in Figure 8. Net exposure times for the high and low-flux spectra for each XIS (and 
the PIN) were approximately 62.1 (46.3) and 71.3 (59.5) ks, respectively. The average 2-10 
keV fluxes were 2.47 and 2.16 xlO -11 erg cm -2 s -1 , respectively. We applied the best-fitting 
baseline model from the time-average spectrum to both spectra. All narrow Gaussian energy 
centroids and widths, (3, i and R- m for the disklines, and log(f) for the high-ionization absorber 
were kept fixed at their time-averaged values. However, all freed parameters were consistent 
at the 90% confidence level. For instance, we find no strong evidence that the column density 
of the high-ionization absorber varies on short timescales. A broadband observation of NGC 
3516 spanning a larger flux range is thus needed to potentially distinguish determine if the two 
power-laws vary in concert. 

6. The Search for red and blue-shifted narrow lines 

We searched for additional narrow absorption or emission features in the Fe K bandpass, 
as seen so far in several Seyferts, including NGC 3516. For instance, Turner et al. (2002) found 
emission lines near 5.57, 6.22, 6.41, 6.53 and 6.9' keV in the November 2001 Chandra-KETGS 
and XMM-Newton EPIC spectra of NGC 3516. One interpretation was that they were red- 
and blue-shifted Fe K lines associated with transient “hot spot” emission on the inner accretion 
disk. 

We searched for such features by adding a Gaussian component to the time-averaged 
spectrum, with width fixed at 0.5 eV, sliding it over 4-9 keV in energy. In addition to the 
aforementioned absorption lines at 6.70 and 6.96 keV, there was only one “candidate” feature 
with Ax 2 < -5.0, an absorption line near 6.0 keV. However, we performed Monte Carlo simu- 
lations to assess the statistical significance of this feature (see §3.3 of Porquet et al. 2004 and 
§4.3.3 of Markowitz et al. 2006 for a description of these simulations), and we found that the 
line was consistent with photon noise. Furthermore, the feature was evident only in one XIS 

2 Indicates that the measured variance is smaller than that expected solely from measurement noise 
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camera, and thus is likely not real. 

As an aside, we note that a 7.47 keV Ni Ka emission line was not detected in the 
time-average spectrum; adding a Gaussian at this energy, we found an upper limit of 3 eV. 

We also extracted time-resolved spectral slices by dividing the time-averaged spectrum 
into five slices 48 ks in duration, with each slice having an exposure time near 27 ks per XIS and 
21 ks for the PIN. We did not investigate longer time slices since they might miss short-lived hot- 
spot emission lines; shorter time slices would have yielded poorer photon statistics. Applying a 
sliding narrow Gaussian over 4-9 keV in each spectral slice revealed only 3 “candidate” emission 
or absorption lines (with Ay 2 < -6). Again, however, the features were seen in only one or 
two XIS cameras, and Monte Carlo simulations showed that not a single candidate feature 
was inconsistent with photon noise at greater than 80% confidence. Analysis of the high-and 
low-flux spectra similarly yielded no significant narrow emission or absorption lines (even at 
6.70 and 6.96 keV). Typical upper limits to the eV of an emission line at 5.57 keV (one of 
the energies of the transient lines in Turner et al. 2002), were <10-15 eV in the time-averaged 
spectrum or any of the sub-spectra. 

7. Discussion 

During the late 1990’s, NGC 3516 typically displayed a 2-10 keV flux of ~4-6 x 10~ n 
erg cm -2 s -1 (e.g., Markowitz, Edelson 2004). During the 2001 XMM-Newton/ Chandra obser- 
vations, however, the observed 2-10 keV flux was much lower: 1.6-2. 3 x 10~ u erg cm -2 s -1 
(Turner et al. 2005). Table 4 lists the inferred absorption-corrected 2-10 keV nuclear fluxes from 
the XMM-Newton observations, as well as during the 2005 Suzaku observation. The Suzaku 
observation apparently caught the source in a similar low level of nuclear flux as the 2001 ob- 
servations. The observed 0.5-2.0 keV flux during the Suzaku observation, however, was ~2-3 
times lower than during the XMM-Newton observations, indicating that the source was still 
heavily .obscured, and confirming that the complex absorption in this source cannot be ignored 
when fitting the broadband spectrum and modeling diskline components. 

7.1. Power-law Components 

The primary power-law observed in the hard X-rays is likely emission from a hot corona 
very close to the supermassive black hole, as seen in all Seyferts. The nature of the soft power- 
law component, however, is not as clear. It could represent nuclear emission scattered off 
optically-thin material (e.g., in the optical Narrow Line Region). The normalization of the soft 
power-law relative to that of primary power-law is 4.2±0.4%, indicating a column density for the 
scattering material of roughly 5xl0 22 cm -2 (assuming a covering fraction of unity), consistent 
with this notion. Using Chandra-ACIS, George et al. (2002) found the extended circumnuclear 
gas to have a 0. 4-2.0 keV flux of roughly 10—14 erg cm -2 s -1 . However, that emission was 
studied over an annular extraction region 3" to 10"; this flux value is likely a lower limit to 
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the 0. 4-2.0 keV flux that the XIS has observed, and so we cannot rule out the hypothesis that 
the soft power-law component is scattered emission. Alternatively, the soft power-law could be 
unobscured, “leaked” nuclear emission as part of a partial covering scenario. In this case, the 
primary absorber would obscure 96% of the sky as seen from the nuclear continuum source. The 
very low level of 0. 3-1.0 keV variability observed could argue for the soft power-law to originate 
in scattered emission, since we might expect to observed variability of the same amplitude as 
the 2-10 keV band only if the soft power-law were leaked nuclear emission. However, this is far 
from certain, as the presence of the soft emission lines in the XIS spectrum could contribute to 
dilution of the observed variability of the soft power-law. A broadband observation spanning 
a larger observed flux range is needed to clarify this issue. The soft power-law could of course 
represent a blend of scattered emission plus leaked nuclear emission. We therefore conclude 
that primary absorber has a covering fraction between 96-100%. 

1.2. Complex absorption 

We detect two zones of absorption: the primary absorber, with a covering fraction of 
96-100%, and the high-ionization absorber, assumed here to have a covering fraction of unity. 
The high-ionization absorber is likely the same as that reported by Turner et al. (2005); we 
fin d a column density TVh of 4.0tf;® xlO 22 cm -2 , consistent with the column density of 2xl0 22 
cm -2 used by Turner et al. (2005). Previous studies of NGC 3516, such as Netzer et al. (2002), 
have discussed in detail the UV absorber, responsible for H Lya, C IV and N V absorption 
features in Hubble Space Telescope spectra (Kraemer et al. 2002). In the X-ray band, discrete 
features associated with Mg VII-IX and Si VII-IX are expected from this component, but 
with the CCD resolution and the calibration-related artifacts near 1.7-1. 8 keV in the XIS, such 
features would likely not be detected by Suzaku. 

Suzaku has found the primary absorber of the hard X-ray continuum to be lowly-ionized 
(log(£) = 0.3±0.1 erg cm s -1 ), with a column density Afl of 5.5±0.2 xlO 22 cm -2 . It is possible 
that it is the same absorber that Turner et al. (2005) designated as the “heavy” partial-covering 
absorber. In this case, the column density has decreased by a factor of 4.5 over a four year 
timescale, while the covering fraction has increased from ~40-60% to 96-100%. However, 
because we have not actually observed a specific blob of absorbing gas moving into or out of 
the line of sight, it is not clear whether the covering fractions derived are associated with single, 
large blobs partially blocking the line of sight to the X-ray continuum source, or if the absorber 
consists of numerous, discrete blobs. On the other hand, given the four year gap between 
the XMM-Newton and Suzaku observations, it is certainly plausible that the 2001 and 2005 
observations could have caught two independent, discrete blobs with differing column densities 
and differing physical sizes and/or radial distances lying on the line of sight. 

To estimate the distance r between the central black hole and the absorbing gas, we can 
use £ = I/2-2oo/(^ r2 )> where n is the number density. L 2 - 200 is the 2-200 keV illuminating 
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continuum luminosity, and £ is measured to be 2 erg cm s _1 . We estimate the maximum possible 
distance to the material by assuming that the thickness A r must be less than the distance r. 
The column density Nn = nAr, yielding the upper limit r < T 2 - 200 / (-/%£)• We estimate the 
2-200 keV flux from the baseline model to be 1.47xlO~ 10 erg cm -2 s _1 , which corresponds to 
L 2-200 = 2.6 x 10 43 erg s _1 (assuming H 0 = 70 km s _1 Mpc -1 and A 0 = 0.73). r is thus <2 
xlO 20 cm (250 light-years), a very loose upper limit encompassing both distances associated 
with the BLR (~10 light-days; Peterson et al. 2004) and a possible cold molecular torus at a 1 
pc radius. 


7.3. Fe K emission components and Compton reflection 


A consequence of Suzaku ’s broad bandpass is that we have deconvolved the broadband 
emitting components, and determined that 1) the broad Fe line was robust in that it was 
required in all models for an adequate fit, and 2) a partial covering component could not 
mimic the curvature associated with a relativistic broad line. The best-fit disk inclination was 
typically <25°. The inner radius was typically <5 R g , suggesting that the black hole is not a 
Kerr black hole, and may be spinning. The line energy was seen to be consistent with neutral 
to mildly-ionized Fe (up to Fe ~ XX; Kallman et al. 2004). 

We measured a narrow Fe Ka line width of 2800ii5o° km s -1 FWHM. We note that 
this velocity is roughly -commensurate with the optical Broad Line Region; Peterson et al. 
(2004) reported FWHM velocities for the Ha and H f3 lines of 4770±893 and 3353±310 km s _1 , 
respectively. The line energy was also consistent with neutral Fe. 

It is possible that the same material that absorbs the hard X-rays is responsible for 
producing the narrow Fe line. The fact that no Compton hump or 6.2 keV Compton shoulder 
are seen, implying an origin for the Fe Ka line in Compton-thin material, supports this notion. 
On the other hand, the material cannot have a column substantially less than 10~ 22 cm -2 or 
else there would be insufficient optical depth to produce a prominent Fe K line. As an estimate 
of the Fe Ka equivalent width expected in this case, we can assume an origin in optically-thin 
gas which completely surrounds a single X-ray continuum source and is uniform in column 
density, and use the following equation: 


EW ca ,\c = 


fcUfaccA 


P(E line) 


( 1 ) 


Emission is assumed to be isotropic. Here, / c is the covering fraction, initially assumed 
to be 1.0. to is the fluorescent yield, 0.34 (Kallman et al. 2004). /k q is the fraction of photons 
that go into the Ka line as opposed to the K{3 line; this is 0.89 for Fe I. A is the number 
abundance relative to hydrogen. We assumed solar abundances, using Lodders (2003). P(E) is 
the spectrum of the illuminating continuum at energy E; E\- ine is the Ka emission line energy. 
Oph(-E') is the photo-ionization cross section assuming absorption by K-shell electrons only 
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(Veigele 1973 s ). 

For TVh = 5.5 xlO 22 cm -2 , EW C& \ C = 29 eV, substantially lower than the observed EW. 
We conclude that it is possible for the absorbing material to contribute to the observed EW, 
but there is also likely contribution from some other (non-continuum absorbing) material, likely 
with column densities >10 23 cm -2 (e.g., Matt et al. 2002). For instance, if the putative cold 
molecular torus does not intersect the line of sight, it could contribute to the observed EW. 
However, the fact that we have resolved the narrow line and obtained a lower limit of 1300 km 
s -1 to the FWHM velocity width could indicate that the torus, where the FWHM velocities 
would likely be < 1000 km s _1 , does not dominate in its contribution to the narrow line EW. 

Finally, we discuss the origin of the material that gives rise to the observed strong 
Compton reflection hump. The primary absorber lacks the necessary column density, and is 
excluded. As far as an origin in the same material that yields the broad or narrow Fe lines is 
concerned, we note that George, Fabian (1991) calculated that R= 1 corresponds to an observed 
line EW of 150 eV for neutral Fe. In our baseline model, we found a best-fit value of R= 1.7, 
which corresponds to a line EW of 255 eV, closer to the observed EW of the broad line (287 
eV in the baseline model) than the narrow line (123 eV). This is consistent with the broad line 
originating in a Compton-thick disk. Combined with the lack of a prominent Compton shoulder, 
the narrow line is thus consistent with arising in Compton-thin material commensurate with 
the BLR. If there does exist a molecular torus, it might not be highly Compton-thick or cover 
a large fraction of the sky as seen from the central black hole, otherwise it would also make a 
significant contribution to the observed Compton reflection strength. 

Suzaku has observed other Seyferts to display reduced levels of variability in the PIN band 
compared to the 2-10 keV band, e.g., in MCG-6-30-15 (Miniutti et al. 2007). This behavior is 
thought to be caused by the presence of the relatively non-varying Compton reflection hump, 
which dilutes the observed >10 keV variability of the power-law component. We would need 
to observe NGC 3516 over a larger >12 keV X-ray flux range than in the Suzaku observation 
to potentially observe any significant variability in the PIN band. 

8. Summary of Main Results 

We have reported on a 150 ksec observation of NGC 3516 made with the Suzaku ob- 
servatory in October 2005. The good exposure times after screening were 135 ks for the XIS 
cameras and 106 ks for the HXD-PIN. 

Our best-fit broadband model included a primary power-law with photon index 
F=1.912±0.025 in our baseline model, absorbed by a column of material with Ah = 
5.5±0.2xl0 22 cm -2 and with log(£)=0.3 erg cm s -1 . We also modeled a soft power-law compo- 
nent which could be scattered nuclear emission off optically-thin material, unobscured “leaked” 

3 http://www.pa.uky.edu/~verner/photo.html 
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nuclear emission, or a blend of both. We conclude that the hard X-ray absorber could be a 
partial coverer, with a covering fraction 96-100%. If this absorber is the same “heavy” ab- 
sorption component reported by Turner et al. (2005) in the 2001 XMM-Newton observations, 
then between 2001 and 2005 the column density of this absorber decreased by a factor of 4.5, 
while the covering fraction increased substantially, leading to an observed 0. 5-2.0 keV flux a 
factor of 2-3 lower in 2005 than in 2001. We also modeled a highly-ionized absorber with a 
column density Ah of 4.0jl3;f xlO 22 cm -2 , ionization parameter log(£) = 3.7lo;f erg cm s _1 , 
and assumed to have a covering fraction of unity. 

Our baseline model also included a dozen narrow emission lines originating in He- like N, 
O, Ne and Mg, H-like C, N, 0, Ne and Mg and at least three RRC lines, consistent with an origin 
in photo-ionized material. However, we cannot exclude a contribution from collisionally-ionized 
material, as suggested by the presence of an Fe L line near 0.83 keV. 

The broad line has been robustly detected: thanks to Suzaku’s broad bandpass, we 
can distinguish between the curvature in the observed continuum due to a partial coverer and 
that due to a broad diskline; we conclude that for this observation of NGC 3516, a diskline is 
preferred and that a partial coverer cannot mimic the curvature of the diskline. The broad and 
narrow lines are decoupled thanks to the narrow response of the XIS and the subsequent high 
signal/noise ratio in the narrow line. Given the observed equivalent widths of the broad and 
narrow Fe lines, the value of R = 1.7t°;5 is more consistent with an origin in the same material 
that gives rise to the broad line, e.g., the accretion disk. The narrow Fe line, meanwhile, has 
a FWHM velocity width of 2800iisoo km s -1 , indicating an origin in material with the same 
velocity as NGC 3516’s BLR. • 
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teams and operations staff for making this observation possible and for assistance with data 
calibration and analysis. This research has made use of HEASARC online services, supported 
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Table 1. Soft X-ray Emission Lines 


Line 

Line 

Intensity 

EW 


F- test 

Energy (keV) 

Identification 

(10~® ph cm" 2 8-1) 

(eV) 

-A% 2 

Probability 

0.366+ 

C VI 

88±26 

32±9 

22.71 

2-lxlO- 5 

0.427+ 

NVI 

42±11 

20±5 

25.00 

7.3x10-® 

0.495±0.005 

N VII/C VI RRC blend 

34±7 

22±5 

47.34 

2.2xlO -10 

0.563±0.004 

O VII 

77±10 

63±8 

154.92 

2.9xl0 -31 

0.665t°;°J2 

O VIII 

15±4 

16±4 

30.04 

6.9xl0 -7 

0.739+ 

O VII RRC 

6.3±2.9 

8.6±4.0 

12.62 

5.4xl0~ 4 

0.830±0.017 

Fe L XVII 

6.2t™ 

8-4t3°7 6 

14.36 

l.lxlO -3 

0.871+ 

O VIII RRC 

13±f 0 

20±f 5 

6.94 

l.OxlO -2 

0.917±0.007 

Ne IX 

16 tj 

26^5 

141.41 

l.OxlO -28 

1.020+ 

Ne X 

5.6±1.5 

14±4 

37.45 

2.9xl0 -9 

1.196+ 

Ne IX RRC 

5.2±1.2 

18±4 

48.22 

1.7xl0 -11 

1.351±+0.006 

Mg XI 

6.9±?;| 

30t| 

103.62 

1.8xl0 -21 

1.471+ 

Mg XII 

2.8±1.0 

14±5 

20.25 

1.2xl0 -5 


Results are for our best-fit baseline model, 
t denotes a fixed parameter. 
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Table 2. Best-fit parameters for the 1-absorber model and the baseline (2-absorber) model 




1-absorber Model 

Baseline (2-absorber) Model 

x 2 / do f 


1521.46/1419 

1485.04/1417 

r 


l-867l°;°42 

1.904±0.025 

Primary Power-law Norm. 

l 

GO 

bo 
l + 
o o 
io it*. 
X 
)— 1 
O 
1 

CO 

9.5±0.5 xlO -3 

Primary absorber 

TVh (10 22 cm” 2 ) 

5.5±0.2 

5.5±0.2 


log(0 (erg cm s -1 ) 

0.3±0.1 

0.3±0.1 

Soft Power-law Norm. 1 


3.9±0.2 xlO -4 

4.0±0.2 xlO -4 

High-ionization absorber 

N& (10 22 cm*” 2 ) 

- 

4-0 Hi 


log(£) ( er S cm s" 1 ) 

- 

o y+0.3 
o>/ -0.7 

Fe Ka diskline 

Energy (keV) 

6 40 +0 - 18 

o.w _ 0 qo 

o 40 +O.I 3 

O.3:0_o.08 


Rin (Rg) 

<3.5 

<3.5 


Inclination i 

<23° 

25 ±1° 


Intensity (ph cm" 2 s" -1 ) 

i. 33+8;55 xio" 4 

1.22lo.io xlO -4 


EW (eV) 

296t|e 

287tH 

Narrow Fe Ka line 

Intensity (ph cm" 2 s” 1 ) 

5.3±0.4xl0 -5 

5.2±0.3xl0 _s 


EW (eV) 

132±10 

123±7 

Compton reflection strength R 

P+-0.3 

1 y+0.4 
L '-0.5 

1 Units of power-law normalization are ph keV -1 cm -2 

s* -1 keV"" 1 at 1 keV. 



2 PEXRAV inclination angle set at 20°. 
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Table 3. Best-fit parameters for models with partial covering components 



Model PCI 

Model PC2 

Model PC3 

X 2 /dof 

1485.69/1416 

1480.48/1414 

1676.57/1418 

r 

1 Q21 +0.027 
i.yoi_ 0 057 

i 079+0.011 
i * y * z -0.036 

9 0^1 +0.046 

4-*OOi_Q 04g 

Primary Power-law Norm. 1 

7.5±0.2 xlO" 3 

9.8±g;^ xlO" 3 

i-3±g;| xio- 2 

Primary absorber iV** (10 22 cm” 2 ) 

5.5±0.1 

6-lio.i 

5.3±0.2 

log(f) (erg cm s' 1 ) 

0.3±0.1 

0.3±0.1 

o.i±g;f 

Partial covering Component IVh (10 22 cm” 2 ) 

5.5 2 

i-o±o:i 

7.8i? 8 s * 

log(£) (erg cm s" 1 ) 

0.3 2 


-0.3±°; 8 7t 

Power-law Norm. 1 

2.4±0.2 xlO -3 

xlO -4 

<2.4x lO" 3 

Soft Power-law Norm. 1 

4.0±0.1 xlO" 4 

3 -7±S:l xl0_4 

4.1±g;S xlO- 4 

High-ionization absorber IVh (10 22 cm*” 2 ) 

4.1±2.5 

3-3±i;| 


log(£) (erg cm s _1 ) 

3.3 ±g;f 

o r+0.8 
°-°-0.2 

3.2±0.1 

Fe Ka diskline Energy (keV) 

6.45±°:JL 

6- 45 t° 0 ; 1 o°s 


N'm (Rg) 

<3.2 

<4.3 


Inclination i 

26ti° 

26±|° 

- 

Intensity (ph cm -2 s" 1 ) 

1.18±g;ff xlO" 4 

1.13±g;f| x 10- 4 

~ 

EW (eV) 

268lg 

2661?? 

- 

Narrow Fe Ka line Intensity (ph cm -2 s -1 ) 

5-4±g;| xl0~ 5 

5-4±g;| xlO" 5 

8.8±0.7 xlO” 5 

EW (eV) 

153±it 

134 t\l 

226±18 

Compton reflection strength R 

-j 7+ 0.2 

i-s±8;l 

2 4 +0-93 

An asterisk (*) denotes that the parameter uncertainty pegged at the hard limit. 


1 Units of power-law normalization are ph keV” 1 cm” 2 

2 parameters tied to those of the primary absorber. 

3 PEXRAV inclination angle set at 25°. 

s 1 keV 1 at 1 keV. 



Table 4. Comparison to the 2001 XMM-Newton observations 



Suzaku 

XMM-Newton 

XMM-Newton 


Oct. 2005 

Apr. 2001 

Nov. 2001 

Absorption-corrected 2-10 keV flux (erg cm” 2 s” 1 ) 

3.4X10" 11 

2.7xl0- u 

1.9xl0- n 

Observed 0. 5-2.0 keV flux (erg cm” 2 s” 1 ) 

1.3x10“ 12 

4.3X10" 12 

2.9xl0- 12 

Primary absorber IVh (10 22 cm” 2 ) 

5.5±0.2 

25±1 

25±1 

Primary absorber covering fraction 

96-100% 

44±6% 

58±5% 

High-ionization Afe (10 22 cm” 2 ) 


1.6+ 

1.6+ 


Results for the XMM~Newton spectra were taken from Turner et al. (2005). 
Results for Suzaku are from the best-fit baseline model. 

T denotes a fixed parameter. 
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10 20 50 


Observed Energy (keV) 

Fig. 1. HXD-PIN spectrum. The upper panel shows the net source spectrum ( gray points ), the back- 
ground ( lower black points ), and the total (source + background) spectrum ( upper black points). The PIN 
spectra have been binned such that the net spectrum has a minimum signal-to-noise ratio of 8cr per bin. 
The lower panel shows the ratio of the net source spectrum to the total spectrum. 
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1 10 

Energy (keV) 


Fig. 2. Data/model residuals for selected broadband models. Panel a) shows residuals to a simple 
power-law. Panel b) shows the residuals after the primary absorber and narrow Fe lines are included 
in the model. Panel c) shows residuals after the soft power-law is included. In panel d), the soft emis- 
sion lines have been modeled. In panel e), a PEXRAV reflection component has been added. In panel f), 
diskline components have been added. In panel g), the high-ionization absorber has been added to yield 
our baseline model. Black points <12 keV denote the (co-added) XIS-FI spectrum. Gray points denote 
the BIS-BI spectrum. Black points >12 keV denote the PIN data. Rest-frame energies are shown. All 
data have been rebinned with a binning factor of 5, though XIS data >7.2 keV have been rebinned by a 
factor of 25 for clarity. 
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Energy (keV) 


3. Ratio of the soft band data to a simple power-law, showing th 
;y, only the XIS BI spectrum is shown. Rest-frame energies are sho^ 


Ne IX RRC 





4 6 8 

Energy (keV) 


Fig. 4. Data/model residuals to various models are shown, focusing on the Fe K bandpass. Panel a) 
shows the residuals when the primary absorber is the only zone of absorption modeled. Panel b) shows 
the residuals to the baseline model, with two zones of absorption; note that residuals near 6.9 keV are now 
improved. In panel c), the broad and narrow Fe lines have been removed from the baseline model. Panel 
d) illustrates that when the broad lines are removed from the baseline model and the model is re-fit, large 
data/model residuals appear in the Fe K bandpass. 
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2.5 



Fig. 5. Left : Contour plot of the Compton reflection fraction R versus the photon index of the primary 
power-law. The black contours show results when the PIN/XIS-FI normalization was kept fixed at 1.13. 
The gray contours show results when the PIN/XIS-FI normalization was left free. Right : Contour plot 
of R versus the PIN/XIS-FI normalization; T was a free parameter. In both plots, dotted, dot-dashed, 
dashed, and solid lines denote 68, 95.4, 99.73%, and 99.99% confidence levels, respectively. 
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Fig. 6. Unfolded broadband spectrum for the baseline model. The thick black line is the total spectrum. 
Thin black lines denote the primary, absorbed power-law, the soft power-law, and the reflection component 
peaking at 20 keV. Gray dashed lines denote the narrow Fe lines. Gray solid lines in the Fe K bandpass 
are the diskline components, gray solid lines in the soft band denote the soft emission lines. 
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Fig. 8. Orbitally binned XIS and HXD light curves. The top and middle panels respectively show the 2-10 
keV and 0. 3-1.0 keV count rate light curves, summed over all four XIS cameras and background subtracted. 
The bottom panel shows the 12-76 keV PIN light curve, also or bitally-binned and background-subtracted. 
The error bars show only the statistical error of the count rates; systematic uncertaintly associated with 
subtraction of the non-X-ray background in the PIN is roughly 0.3 ct s -1 . The mean count rates are 2.41 
ct s _1 for 2-10 keV, 0.20 ct s -1 for 0. 3-1.0 keV, and 0.16 ct s -1 for 12-76 keV. 
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